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Abstract. We consider the radiation field operators in a cavity with varying dielectric medium in 
terms of solutions of Heisenberg's equations of motion for the most general one-dimensional quadratic 
Hamiltonian. Explicit solutions of these equations are obtained and applications to the radiation 
field quantization, including randomly varying media, are briefly discussed. 
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1. Canonical Quantization 



Radiation field quantization in the vacuum was introduced in original works of Born, Heisenberg 
and Jordan [ID] , [TT] (see also books on quantum electrodynamics [I], [3J, [6] and quantum optics 



6], [63], [82], [HI], [93J ) . A modern mathematical approach to quantization of mechanical systems 
is discussed in detail, for example, in [9], [26], [85] . and/or [90] (see also [45] and the references 
therein). For a classical Hamiltonian system one replaces canonically conjugate coordinates and 
momenta by time-dependent operators q\{t) and P\{t) that satisfy the commutation rules 

[q x (t), qfl (t)} = \p x {t), pM = 0, [q x (t), p„(t)) = ihSxp. (1-1) 

00 ■ 

{Sj ■ The time-evolution is determined by the Heisenberg equations of motion [35j : 

R; j t Px(t)= t -[p x (t),'H], f t qx(t) = ^[q x (t),n], (1.2) 

O ■ with appropriate initial conditions^! 

Traditionally, the electromagnetic-field quantization is considered under the assumption that the 
field occupies an empty box PQ, [3], [20], [27]. The quantization of the field in a uniform dielectric 
medium was discussed in [25], [37], [35], [10], [2], [UJ (see also [2D], [21], [53] and the references 
therein). Yet the problem of electromagnetic-field quantization in time-dependent nonuniform linear 
nondispersive media remains an active research topic up to now [2], [7J, [12], (T3J, [20] , [30], [36] . 
[38], [H3. EFJ. [M]. |7J. [75], [76], [77], [88], [52]. 

In the present letter, we study the radiation field operators in a varying dispersive medium, 
which is mathematically described by the most general phenomenological quadratic Hamiltonian 
H in an abstract Hilbert space. (We concentrate on a single photon cavity mode, say v, with 
frequency u v = 1 and use the units c = h = 1. From now on, we shall usually omit the indices when 
dealing with the single mode under consideration.) In particular, our approach gives a natural 
description of squeezed photons that can be created as a result of parametric amplification of 
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l The standard form of Heisenberg's equations can be obtained by the time reversal t — > — t. 
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quantum fluctuations in the dynamic Casimir effect jTTJ, [18], [51], [67], [69], [83], [M] and/or by 
similar dynamical amplification mechanisms including the Unruh effect [91J and Hawking radiation 
[8], [33], [31]. It is also useful for quantum fields propagating in nonstationary external potentials 
[7J, [60], [92], [83], and for photon quantization in randomly varying media. 



2. Solution of Heisenberg Equations for Nonautonomous Quadratic Systems 

Our main result is the following. 

Theorem 1. The solution of the Heisenberg equations of motion ( II. 2p for the nonautonomous 
quadratic Hamiltonian 

H = a{t)p 2 + b{t)x 2 + c{t)xp - id{t) - f{t)x - g{t)p (2.1) 
(a, 6, c, d, f, and g are suitable real-valued functions of time only) has the form 

Pi(t) = Mz|W, 9A(t) = M±fW. ( , 2) 

i/ere, i/ie time- dependent annihilation b(t) and creation b\t) operators are given by the Ansatz 

T . . e -2n(*) / p _ 2a (t) x - 8 (t) 

m = ~yr r (t) 1 + £ (i) + ' — W) 

r+/ x e 2 ^W/^/x .p - 2a (t) £ - o" (t) \ 

&t(t) = — ^- (t) x + e (t) - ^ ^ Lij (2.3) 

m terms of solutions of the Ermakov-type system 

^ + 6 + 2ca + 4aa 2 = a/3 4 , (2.4) 
d/3 

+ (c + 4aa) /3 = 0, (2.5) 

^ + a/3 2 = (2.6) 

anc? 

dS 

— + {c + 4aa)5 - f -2ga = 2af3 3 e, (2.7) 

Cot- 

- - (0 _ 2*5) = 0, (2.8) 

^-^o + ao 2 = a/3V. (2.9) 

T/ie time-independent (self-adjoint) operators x andp obey the canonical commutation rule [x,p] = i 
in an abstract (complex) Hilbert space which implies that the relation 

b(tft(t)-ft(t)b(t) = l (2.10) 

holds at all times. 
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Proof. These results can be verified by a direct, but somewhat tedious, calculation when one expands 
the solution in generators of the Heisenberg- Weyl algebra, namely {l,x,p}, with undetermined 
time-dependent complex coefficients and simplifies the commutators. The substitutions f l2.2j) -( l2~3l) 
allow us to derive equations (I2.4l) -f [2~8l) . say from the first Heisenberg equation (12. lip below. (A 
Mathematica based proof is available on the article's website; see notebook HeisenbergOscillators.nb 
and jl3] for an important program ingredient @) Equation (12.91) . which determines the global phase 
of the corresponding Fock states in the Schrodinger representation, does not show up in this proof, 
but will appear later (see Lemma 2). □ 

By back substitution, we see that the operators b(t) and o'(t) are solutions of the Heisenberg 
equations 

d b( t ) = % \b(t), H] , H\t) = % \V(t), H} , (2.11) 



dt 

subject to the initial conditions 



dt 



-2vr(o) / .p - 2a(0)x - 6(0) 

2*7(0)/ _ _ 2a(0)x - S(0)\ 

(0) = ~7T \ m)x + £{0) ~ W) — J • ( } 

To a certain extent, the creation and annihilation operators (12. 3p allow us to incorporate the 
Schrodinger symmetry group of the harmonic oscillator, originally found in coordinate representa- 
tion [70], [71], into a more abstract Heisenberg picture. (For the sake of simplicity, we have restricted 
ourselves to a single photon mode v with frequency oj v = 1; see [49J for a detailed investigation of 
the special case of uniform media.) 

A concept of dynamical invariants for generalized harmonic oscillators, which is crucial for con- 
structing the corresponding Fock states from our creation and annihilation operators, has been 
recently revisited in [T5], [SI], and [BE] (see [2TJ, [23], [SS], [H5], [66J and the references therein for 
classical accounts). 

3. Solving The Ermakov-Type System 

A general solution of H2.4)) — f!2.9[) is given by Lemma 3 of [51] in a real form (see also [H] and [61]). 
In order to proceed to a more compact form, one needs to recall some notation. The substitution 

1 p! d 

a = (3.1) 

4a n 2a 

reduces the inhomogeneous equation (12.41) to the second order ordinary differential equation 

// - r(ty + Aa(t)fi = co(2a) 2 /?V, (3-2) 
which has the familiar time- varying coefficients 

T ( t ) = --2c + 4d, a(t) = ab-cd + d 2 + -(--^] . (3.3) 
a 2 \ a a J 

(In ( 13 .2p and in the rest of the paper, we use a formal 'binary' parameter cq = 0, 1 for the sake of 
convenience.) 



2 See also Koutschan.nh 
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The time-dependent coefficients a , (3 , 7 , 5 , e , k , which satisfy the homogeneous (Riccati- 
type) system $H$-<&M, are given by (cf. [it], P], [SE]) 

Mt) = 4t#| - (3-4) 
ow 4o(*)/i (t) 2a(t)' v ' 

f3 Q (t) = -^Pr, \{t) = exp (- [ (c(s) - 2d(s)) rfs^) , (3.5) 



/i (t)' 



7o(t) - 2/ , 1 (0)Mt) + 2a(0)' (3 - 6) 



and 



M*) Jo l\ a(s) J 2a{s) 



ds 



\{s) 



(3.7) 



2a(t)X(t) . . . „ a(sWs)A(s) , . 
/w) (M s )) 



MO io (M*))' 

-2 t4fs Ms)So(s)) (f(s) - ^f\g(s)) ds (3.9) 
7o M s ) V a(s) / 

(5o(0) = — £o(0) = S"(0)/ (2a(0)) and Kq(0) = 0), provided that /i and //j are the standard (real- 
valued) solutions of equation (13. 2p when Co = corresponding to the initial conditions /i (0) = 0, 
/i' (0) = 2a(0) 7^ and /MO) 7^ 0, /^(O) = 0. (Proofs of these facts are outlined in [H] and [T6] . 
The integrals are treated in the most general way, which may include stochastic calculus; see, for 
example, [28], [73] and the references therein.) 

Here, we would like to present a new compact form of these solutions. Let us introduce the 
complex- valued function 

z(t) = (2a(0) + ^) M0 + ^ + */3 2 (0K(t) (3.10) 

V a (°)/ MO) 

(a complex parametrization of Green's function and linear invariants of generalized harmonic oscil- 
lators are also discussed in [22], [23] and [31]). Then 

z(t) = c x E{t) + c 2 E*(t), (3.11) 

where the complex-valued solutions are given by 

E(t) = ^ + ilh®, E*(t) = ^ - z> (t), (3.12) 

M0) MO) 

and the corresponding complex-valued parameters are defined by 



ci =^^- ! ( a(0)+ ^)J- c2= ^^ +! r 0)+ ^)J' (3 ' 13) 

with 

Cl + c 2 = l, |ci| 2 - |c 2 | 2 = ci -c* = /3 2 (0). (3.14) 
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In addition, 

z(0) = c 1 + c 2 = l, z'(0) = 2ia(0) (ci - c 2 ) . (3.15) 

When written in terms of the complex function z in (I3.10p . the complex conjugate functions E and 
E* defined in (I3.12p become 

Cl - C 2 Ci - c 2 



and 



z-z* (cl -c* 2 )z + (ci - c 2 ) Z* 

^° 2<(c 1 -c5)' ^(0) 2( Cl -c*) • 1 ' ] 



One can readily verify that 



1 (z-z*)' d 



"o = - A — + 7r> 0o = — 2aA 



Ci — c 2 



4a z-z* 2a z — z* 

(Cl-^)2 + ( Cl -C 2 )^ d(0) 

70 2. (z-z*) + 2aT0)' (3 - 18) 

and equations (I3.7p -( 13T9~1) can also be rewritten in terms of the function z in view of (13.171) . Finally, 
we introduce a second complex function, 

C{t) = e(O)0(O) + i (5(0) + e (t)) =c 3 + te , c 3 = e(0)/3(0) + z<f(0), (3.19) 

and indicate the inverse relations between the essential, real and complex, parameters: 

«(0) = ^ " ^y, /3 2 (0) = Cl - cS = | Cl | 2 - |c 2 | 2 , (3.20) 

and 

6(0) = *=A, £ (0) = ± C3 + C 3 . (3.21) 

Zi% c\ I \ 1 2 i 1 2 

2y |ci| - |c 2 | 

Then the solution of the initial value problem for the Ermakov-type system can be expressed in 
terms of the complex function z in (13.101) as given in the lemma below. 



Lemma 1. The system (12.41) - (12.91) is solved by 



|2 _ | |2 

« = «o + A 2 ^|^4, /3 = ±A V ' C , . m , 7 = 7(0)-^argz, (3.22) 
2i\z\ z-z* \z\ 2 

J = *o + A C *~ C * 2 ** , e = ± CZ ; +CV , (3.23) 

2llzl 2N v /|c 1 | 2 -|c 2 | 2 

(C^+ , N 

« = «(0) + «o + 7^ i,\ , 2 (3-24) 

8« (ci — Cg) |z| 

wi/i z and C as in (13. lQf) and (13. 19ft . respectively. (The solution of the Ermakov-type equation (13.21) 
is (jwen by n = /x(0) 

Proof. This amounts to a straightforward calculation using Lemma 3 in [51]. □ 
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As a consequence, one gets 

2i (a — «o) z + z* 



/3 2 z- z* 



and 



(3.25) 



i (a - a ) + ^- = /3 2 — Z -—, i (a - a ) - %r = /3 2 - (3.26) 

2 z — z* 2 z — z* 

, .8 -So (z . S-6 Cz* . Q _, 

£+J T = i)R' £ - l HT = W)¥\' ( } 

^-(°)+-o + ^s-^Ms(0). (3.29) 

These "quasi-invariants" can be useful, for example, when making a comparison of calculations 
done by different approximation methods. 

Examples of explicitly integrable quadratic systems are discussed in [15], [16], [22], [23], [56], [57] . 
and [HS] (see also the references therein). 



4. Single Mode Fock States for Nonautonomous Quadratic Hamiltonians 
The time-dependent quadratic operator (see [81] ) 

= - [a(t)a f (t) + a j {t)a{t)] , (4.1) 
with the defining property 

+ EH - HE = 0, (4.2) 

extends the standard Hamiltonian (and/or number operator) for any given solution of the Ermakov- 
type system. (The initial data play a role of integrals of motion and/or quantum numbers for the 
creation and annihilation operators in the Heisenberg representation.) We use the substitution 
b(t) = e~ 2i7 a(t) and 6 f (t) = a){t)e 2i ^ with \a(t) ,a) (t)] = 1. The oscillator-type spectrum, 

E(t)\* n (t))=(n+^\ |¥„(t)>, (4.3) 

can be obtained in a standard fashion (with the aid of modified variable creation and annihilation 
operators; cf. [1], [9]): 

a(t) \* n (t)) = \V n -i(t)) , a\t) \* n (t)) = v^TT \^n+i{t)) ■ (4.4) 

Here and in what follows, it is convenient to use the orthogonality relation (\P m (£), ^ n (t)) = 5 mn A -1 
with /3(0)//(0) = 1. 

Now we can analyze abstract Fock states in the Schrodinger representation. 



E(t) 



(p — 2ax — S)' 
J 2 



+ (fix + e) 
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Lemma 2. Leta(t) \^> Q (t)) = 0- The dynamic Fock states given by 



l^nW)=e l(2ri+1)7+iK |^n(t)) 

satisfy the time- dependent Schrddinger equation 

d 



3 i(2n+l)7+ire 



(at(t))> (t)) 



i- = H |V „(«)> 



(4.5) 



(4.6) 



t/ie general quadratic Hamiltonian (12.11) provided that equation (12.91) /or t/ie global phase holds 
and (<$> n ,d<$> n /dt) = for $ n = A 1/2 e -< aa;2+5x )^ n . 



Proof. From (I4.2p - (j4.3j) . one gets, formally, 



n 



*^ n (t)>-# |^ n (t)> 



(4.7) 



Therefore 



i^\Mt))-H\1> n (t)) = c n \1> n (t)) 



where Imc„ = in view of the normalization condition; see 
vacuum state is nondegenerate.) 

Here. 



(4.8) 

(We assume also that the 



H = ap 2 + bx 2 + - (px + xp) + - (c — 2d) — fx — gp, 
and the position and linear momentum operators are given by 



1 

x = — 



—= (a + a?) — e 



V 



IV2 

,^ v2a ,^ ^ n . 2ae 



In terms of creation and annihilation operators, the Hamiltonian takes the form 



H 



a ( Aa 2 n ~\ b + 2ca 
P 2 + 



2(3 Z 



+ 



a /4a 2 



-/3 2 + 



- (c + Aaa) 

b + 2ca i . . 

-5 V - (c + Aaa) 

2(3 2 2 v ' 



(at) : 



1 

+ 2 
+ V2 



a [(3 2 + 



b + 2ca 



(aa) + a^a) + - (c — 2c?) 



4aa + c 
^X 



2ae 



— 2 (6 + ca) 



4aa + c 
^X 



2ae 



2 (b + ca) 



f + 2ga 

2(3 



f + 2ga 

2(3 



- ^{P{9- 2a5) + e (c + 4aa)) 



(4.9) 

(4.10) 
(4.11) 
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+ a[ 5- 



2ae 







6- 



2ae 



ce 



(4.12) 



where we have corrected a typo in [ST]. As a result, by (|4.4p — (I4.3p we have 



\Re(V n ,HV n ) = in + 



a \(5 2 + 



4a 5 



+ a[5 



+ 



2ae 



b + 2ca 







5 



2ae 

T 



9 + j) (1-13) 



in terms of solutions of the Ermakov-type system. 

In order to complete the proof, one can repeat the evaluation of Berry's phase [I] for generalized 
harmonic oscillators, given in [81] in coordinate representation, in a more abstract form. Indeed, 



(2 + 1)^ + ^ K 
dt dt 



+ (^ n ,i^ n ) = ,). LI ! 



In general, 
and 



A(* r 



dt 



; (l/2)/(c-2d) dt e i(ax 2 +Sx)^ 



/$ $ ) = s 



(4.15) 



i ( $, 



da 
~dt 



d<K 
dt 



dt 1 dt 

where {& n ,d$ n /dt) = by our hypothesis. Here 



l^X (tf n , X 2 ^ n ) + i^-X (9 n , X^ n ) + I (C - 2d) 



As a result, 



A (tt n , x^ n ) = -e/T 1 , A (tt n , x 2 ^ n ) = (3~ 2 (e 2 + n + ~J . 

1\ da x d5 



A Re i ( * 



d^ r , 

dt 



e" + n + 



2 dt 



+ 80' 



dt 



(4.16) 

(4.17) 
(4.18) 



Finally, in view of (|4.8p . (14. 13j) and ( I4.18p . we obtain (after some simplification) that 

c n = ($ n ,id<$> n /dt) = 
for any given solution of the Ermakov-type system (12.4p - (l2.9p . This completes the proof. 



□ 



Remark 1 . In coordinate representation, when $ n is, essentially, the real- valued stationary orthonor- 
mal wave function for the simple harmonic oscillator with respect to the new variable £ = (3x + e 
(see [52], [53], [72], and [81] for more details), the equation ($ n , d$ n /dt) = is valid due to the 
normalization condition ||$ n || 2 = 1- In general, one gets c n = (<& n ,id<& n /dt), and the previous con- 
nection is associated with a transport law for line bundles in the Hilbert space, namely, the change 
d$ n is orthogonal to $ r 



(4.19) 



The last lemma can be reformulated in terms of an analog of Berry's phase [3], 



d6 r . 



— — = A Re ( 7jj w i—^p 



dt 



d 



dt 



d<Pn 

dt 



where <p n (t) = — (2n + l)j{t) is the dynamical phase [81] (see also [87] for an example). 
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Lemma 3. The Fock states, given by (14. 5p in terms of solutions of the Ermakov-type system, satisfy 
the Schrddinger equation ( 14. 6 p if and only if the derivative of the Berry phase is evaluated by two 
equivalent expressions (40) and (48) in [HI] . 



5. Expectation Values and Variances 



By (I4TT0D - (ICTD and one gets 

<*.(«),**.(«)> = —. *=$ = tHt = -!$. <") 

= A (*-^Y ^M^TTTT"^')- 5 ^ < 5 - 2 > 

(the Ehrenfest theorem is discussed in |54j). 
The standard deviations, 

(Sp) 2 = ^1 = W)~ (pf , (Sx) 2 = ^1 = pj _ , (5.3) 



(*p) 2 = ( - + ^ + ^ V (^) 2 =(n + l -)p- 2 (5.4) 



are given by 

(fin) 2 = ( n 4- , 
in terms of solutions of the Ermakov-type system. In particular, 

(Sp) 2 (Sx) 2 = L+^Jfl + ^j >\ (5-5) 

as required by the fundamental Heisenberg uncertainty relation [7J, [15], [35]. The minimum- 
uncertainty squeezed states occur for n = if a (t min ) = 0. 



6. Electromagnetic-Field Quantization in Varying Media 

In the macroscopic approach, one can present the (noncommuting) vector field operators of 
the electric displacement D(r,t) and the magnetic induction B(r,t), which fully describe the 
properties of the quantized electromagnetic radiation inside a cavity filled with linear nonstationary 
dielectric material (with factorized electric permittivity and magnetic permeability tensors [20J ) , by 
the eigenfunction expansions (cf. [7J, [25], [35], [12], [75], [82] ) 

D(r,t) = J2x v (t)p v (t)D v (r), 

V 

B(r,t) = J2^v(t)q v (t)B v (r), (6.1) 

V 

in the Heisenberg picture when the time evolution is introduced through the equations (I2.2l) - (l2.3p . 
which provides a more direct analogy between quantum and classical physics [7J, [32J. For a dis- 
cussion of properties of the stationary orthonormal eigenfunctions D v (r) and B v (r) defined by the 
geometry of the cavity and given boundary conditions, see [2J, [7J, [TJJ, [13], [IS], ED]; ES]) (29J, 

nz], m, m, m, m.- ^ m. 
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In view of the phenomenological Maxwell equations, the single electromagnetic radiation mode 
v in a cavity resonator is analogous to a parametric driven harmonic oscillator [13], [20], [27], [64J, 
[71], [75]. This analogy between classical mechanics and electrodynamics allows one to determine 
functions x v W and ro„ (t) from the electric permittivity, magnetic permeability, and conductivity 
of the (slowly-) varying medium in connection with the quadratic Hamiltonian (12.11) (see Appendix A 
for more details). After quantization of the field Hamiltonian, the time-dependent operators p v (t) 
and q v (t) are determined by Theorem 1, and the corresponding Fock states are constructed in 
Lemma 2. (The average fields obey the classical Maxwell equations [7J.) Methods of stochastic 
calculus |73j can be used in the case of randomly varying media. 

In Schrodinger's picture, for the diagonal matrix elements of the field oscillators, we get 

(D{r,t)) = D v (r) Xv (t) ^ n (t),jnP n (t)) , 

(B(r, t)) = B v (r)w v (t) (Vn(*W n (t)> , (6.2) 

with (I5.1l) - (l5.2p for a single mode v, and the corresponding variances can be obtained with the 
help of (15.41) . In the autonomous case (see [IS]), the variances are given (up to a normalization) by 
equations (A.4)-(A.5) of [52]. 

7. Summary and Applications 

In this letter, the radiation field operators in a cavity with varying dielectric medium are con- 
structed in terms of explicit solutions of Heisenberg's equations of motion. The phenomenological 
quadratic Hamiltonian under consideration corresponds to the most general (single mode) one- 
dimensional mathematical model of quantization in an abstract Hilbert space. Nonstandard solu- 
tions of these equations are obtained and applications to the radiation field quantization, including 
randomly varying media, are briefly discussed. 

For most applications in (nonlinear) optics, the electromagnetic field can be treated classically. 
But, when quantum limits are approached and one is interested in the photon statistics of the field, 
a quantum description is required (see [24], [29], [36], [37J and the references therein). Explicit 
form of the Bogoliubov transformation (12. 3p (for nonautonomous quadratic systems) is one of the 
starting points in this approach [7J, [20], [74], [95]. Interaction of multiple modes and microscopic 
(lossy) medium models are also under consideration [39], [44], [47], [55], [58], [59], [77], [83], [88] . 

The problem of quantization of electromagnetic field in material media remains important in 
view of recent trends in the flourishing cavity QED [IT], [19], [44], [69], [77] and for experiments in 
quantum optics in media [7J, [21], [3D], [23 , [IB], [92J, which may result in a better understanding 
of the interaction of light with matter. Among other possible applications of the electromagnetic 
wave propagation in time-dependent media are the modulation of microwave power [63], wave 
propagation in ionized plasma [SO], and magnetoelastic delay lines [79], [3D] (see also [13] and the 
references therein). 

Acknowledgments. We thank Michael Berry for bringing Ref. [5] to our attention and Luc Vinet 
for a valuable discussion. This research was partially supported by an AFOSR grant FA9550-11-1- 
0220 and by the Austrian Science Foundation FWF, grants Z130-N13 and S50-N13, the latter in 
the framework of the Special Research Program 'Algorithmic and Enumerative Combinatorics" . 
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The phenomenological Maxwell equations in linear, passive, nondispersive, time- varying dielectric 
and magnetic media without sources, namely 

1 BB 

cm\E = — , div£> = 0, (A.l) 

c dt ' ' K J 

cuy\H = - — , div# = 0, (A.2) 

c dt y J 

D = e(r,t)E, B = Jl(r,t)H, (A.3) 



with the help of the vector potential 



B = curlA, E = -~, (A.4) 

c at 



dA\ 

div [e— =0, = 0, (A.5) 



and imposed gauge conditions 



dt J 

can be reduced to the single second-order generalized wave equation 



cml(r 1 cmlA) + ±l(e?£)=0. (A.6) 



c 2 dt\ dt J 

Here, we recall the case of factorized (real- valued) dielectric permittivity and magnetic permeability 
(tensors) , 

e (r, t) = £ (t) e (r) , # (r , t) = 77 (*) £ (r) , (A.7) 

considered in |201. 



The solution of the classical problem for a given single mode, say v, has the form A (r, t) 
u (r) q (t) and 



J B = gcurlu, D = —-— eu (A.8) 



curl I — curl it ) = v 2 eu, div (eu) = 0, 



provided that 



d 2 q £' c 2, u 2 . , 

_ + ? _ + _ g = 0, t, = constant, (A.9) 

and certain required boundary conditions are satisfied on the boundary of the cavity (see [20] . [25] . 
[29]. [4"7]. p] for more details). 

Thus we can choose a = — £, b = —c 2 v 2 / (4£ 2 ?]) and c=d=f=g=0m Theorem 1 for the 
quantization of the mode of the electromagnetic field under consideration. (See also [2], [13], |64j, 
, and [75].) 
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